Anomalous Supercurrent from Majorana States in Topological Insulator Josephson 

Junctions 



Andrew C. Potter and Liang Fu 
Department of Physics, Massachusetts Institute of Technology, Cambridge, MA 02139, USA 

We propose a Josephson junction setup based on a topological insulator (TI) thin film to detect 
Majorana states, which exploits the unique helical and extended nature of the TI surface state. 
When the magnetic flux through the junction is close to an integer number of flux quanta, Majorana 
states, present on both surfaces of the film, give rise to a narrow peak-dip structure in the current- 
phase relation by hybridizing at the edge of the junction. Moreover, the maximal Majorana-state 
contribution to Josephson current takes a (nearly) universal value, approximately equal to the 
supercurrent capacity of a single quantum-channel. These features provide a characteristic signature 
of Majorana states based entirely on supercurrent. 

PACS numbers: 71.10.Pm, 73.23-b, 74.50+r, 74.90.+n 



Majorana bound-states in superconductors are local- 
ized quasi-particles that are equal weight superposition 
of electron and hole, which have non-Abelian braid- 
ing statistics[l, 2}. The presence of Majorana bound- 
states can produce unusual transport phenomena such 
as a 47r-periodic Josephson effect [3, 4], resonant An- 
drccv reflection [5], and inherently nonlocal transport [6] . 
Over the last few years, proposals for realizing Majo- 
rana states in various superconducting solid-state sys- 
tems have sparked tremendous interest and intensive ac- 
tivity (see [7, 8] and references therein). Recent exper- 
iments on semiconductor nanowires with proximity in- 
duced superconductivity[9, 10] observed signatures [11- 
14] interpreted as evidence of Majorana states. 

Among the proposed material systems for realizing 
Majorana states, the superconductor (S)-topological in- 
sulator (TI) hybrid structure[15] has several distinctive 
features. First, as a parent phase for Majorana states, 
the S-TI interface is topologically nontrivial even at zero 
magnetic field. Consequently, the induced superconduc- 
tivity on the topological surface states of a TI is immune 
to disorder[16]. This provides a robust route to realiz- 
ing Majorana states at elevated temperatures. Second, 
by their topological nature, surface states of a TI extend 
throughout the entire sample boundary. The extended 
nature of such states has motivated us to propose a TI- 
based setup for detecting Majorana states. 

In this paper, we study anomalous Josephson cur- 
rent signatures of Majorana states in a superconductor- 
topological insulator thin film-superconductor junction 
under an applied magnetic field. Importantly, the super- 
conductors and the TI film are arranged in such a way 
that induced superconductivity exists on both top and 
bottom surfaces of the TI. The conceptually simplest ex- 
ample of such a setup is shown in Fig. la, in which the 
superconductors wrap around the TI in the entire y and 
z direction and the supercurrent flows in the x direction. 
This geometry will provide a useful starting point for the- 
oretical analysis. However, our results will also apply to 
the simpler to fabricate geometry where the supercon- 




FIG. 1. Panel a) shows a depiction of the conceptually sim- 
plest device geometry considered in the main text (simpler to 
fabricate devices such as those measured in Ref. [17-22] are 
also discussed). The top and bottom surfaces of the topologi- 
cal insulator (grey layer labeled TI) are in contact with super- 
conductors (blue layers). As the global-phase offset between 
the left and right superconductors (6q = 9r — 6l) is adjusted, 
Majorana modes (shown as red-circles) bound to Josephson 
vortices are created at one end of the junction, move along the 
junction, and fuse on the opposite side of the junction. Panel 
b) shows the local phase difference y for along the junction 
3<1>() > "3>b > 2$o and fixed 8o, and the corresponding mass 
term of Eq. 3. Wherever 8 y = ir mod 2n, there is a local gap 
closing that binds a Majorana state. 

ductors are deposited only on the top surface, provided 
that the superconductivity is transmitted to the bottom 
surface through either bulk states (which are present in 
all experimentally realized TI materials), or side surfaces. 
Supercurrent through such SC-TI-SC junctions has been 
recently observed[17-22]. 

Our main findings and the basic physics behind them 
can be stated in simple terms. Applying a magnetic field 
to the S-TI-S junction induces a one-dimensional array 
of Josephson vortices, due to the winding of the phase 
difference between the two superconductors, 9 y cx By, 
in the y direction parallel to the junction. Each vortex 
traps a localized Majorana state on the top surface of the 
TI, and one on the bottom surface. The two are aligned 
vertically (see Fig. 1). The magnetic flux through the 
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junction determines the number of Majorana states and 
their separation. The global phase offset 8 y=a defined 
at a reference point y — is an independent variable, 
which can be controlled in a SQUID geometry. For a 
given magnetic field, increasing 8 shifts the positions of 
all Majorana states along the junction, towards one edge 
of the TI sample. When the first pair of Majorana states, 
one from the top surface and the other from the bottom, 
approach the edge, they become hybridized due to the 
wavefunction overlap on the side surface, and eventually 
annihilate each other. As #o increases further, this fu- 
sion event repeats for the next pair of Majorana states 
and so on. Meanwhile, an opposite event occurs at the 
other edge, where pairs of Majorana states are created. 
Decreasing the phase difference or reversing the direction 
of the magnetic field will reverse the direction of motion 
of Majorana states and interchange the role of the two 
edges. 

Each creation or fusion event causes an energy split- 
ting of the pair of Majorana states involved, and for this 
reason, contributes to the DC Josephson supercurrent 
an amount that is independent of the width of the junc- 
tion. When the magnetic flux through the junction is 
(close to) an integer multiple of flux quanta, the normal 
contribution to the supercurrent oscillates with position 
y along the junction and (nearly) cancels to zero, while 
the Josephson supercurrent Im from splitting Majorana 
states dominates. Advancing 80 by 2ir shifts each Majo- 
rana state to the next one. During this cycle, the first 
pair of Majorana states in the array is fused, and a new 
pair is created and takes the place of the last one. This 
leads to a narrow peak in Josephson current Im as a 
function of phase difference within each 2tt period (see 
Fig.ld). 

Model - Wc now derive the Majorana states and the 
Josephson current-phase relation for a short and wide S- 
TI-S junction, under the condition L < £ < W where 
L is the length of the junction, W is the width, and £ 
is the coherence length of topological surface states with 
proximity-induced superconducting gap. Moreover, for 
simplicity we consider the case that the Josephson cur- 
rent is unable to screen the external magnetic field. 

The Hamiltonian for the top surface of the S-TI-S junc- 
tion is 

[vz ■ (tv x s) — ^i(r)] ip(r) 



+ 



A(rW(r)v4(r) 



(1) 
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(ipf(r), ipi{r)) describes the TI surface states, which have 
a Dirac dispersion with Fermi velocity v. Due to doping 
from the superconductors, the chemical potential in the 
superconducting region fi(r) = [i! at > L/2 is differ- 
ent from the junction area /i(r) = /j, at \x\ < L/2. Aj(r) 
is the vector potential associated with the magnetic field 



B. It is convenient to work in the Landau gauge A x = 0. 
A y (r) is then given by A y (r) = —Bx for |a;| < L/2, 
-BL/2 for x > L/2 and BL/2 for x < -L/2. In this 
gauge, the superconducting gap A(r) is given by 



Ae l0 y/ 2 x > L/2 
A(r) = { Ae~ ie y/ 2 x < -L/2 
\x\<L/2 



(2) 



The phase winding 8 y — Try/Is ensures zero supercur- 
rent along the y direction in the superconducting region. 
Here, the magnetic length is defined as: Is — W^-, 
where $s is the magnetic flux through the junction[24]). 

Since A(r) varies slowly with the position y, we use 
semiclassical method to first solve the Hamiltonian with- 
out the kinetic energy term —ivd y s x at an arbitrary y. 
This one-dimensional problem was previously solved for 
TI Josephson junction[15] and for a related problem in 
graphene[25]. We quote their results below. When the 
junction length L is shorter than the coherence length 
£ = v/A, there is a single pair of subgap Andreev bound 
states at energy ±E(y), where E{y) = Acos{8 y /2) de- 
pends on the local phase difference and oscillates with the 
position y (see Fig. lb). Importantly, the two branches 
of Andreev bound states cross each other at E(y) = 0, 
when the local phase difference 8 y is an odd multiple of 
7r. Such a level crossing corresponds to a local change in 
the fermion parity over a cycle, which is robust against 
perturbations [26] . 

Majorana Bound-States - We now take the kinetic 
energy along the y direction into account, to obtain the 
low-energy spectrum of the junction near E = 0. For 
this purpose, it is instructive to first consider the region 
near a level crossing point y — y\. This region can be 
described[15] by the following effective Hamiltonian in 
terms of two branches of counter-propagating Majorana 
fermions ^l,r- 
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ivMilLdy'iL - iRdy^R) + iE(y)j L ~/ R , (3) 



with E(y) s=s irA(y — yi)/lB, an d where the velocity % 
arises from the kinetic energy along y direction. The 
bound state energy E(y) corresponds to the hybridization 
between the two Majorana states, and vanishes linearly 
at y = 2/1. 

The sign reversal of E(y) as a function of y gives rise 
to a zero-energy Majorana bound state that is spatially 
localized in the junction at the position y . The corre- 
sponding Majorana operator 71 is given by: 



7i 



/ 



dy- 



-(y- yi f/2\ 2 B x 



V2 



bdy) + i R {y)] (4) 



where the decay length in the y direction Xb is: Xb = 
\ZvmIb/ttA. This Majorana state is confined by the TI 
band gap in the z direction, the proximity-induced su- 
perconducting gap in the x direction, and the magnetic- 
field-induccd linear potential E(y) in the y direction. The 
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decay length in the z direction is given by the penetration 
depth of topological surface states into the bulk, which 
is typically a few nanometers. The decay length in the 
x direction is given by the coherence length £ (typically 
a few hundred nanometers). The velocity vm depends 
on intrinsic properties of the junction. For a ballistic 
junction, it was found that vm — v(A//j,) 2 for n = // 
in Ref.[15] and vm — f(A//i) sin(/xL/u) for \i <C jj! in 
Ref.[25]. Since A/it is typically of the order 10 -3 , in both 
cases the decay length of Majorana states X B along the 
junction is much smaller than their separation. Therefore 
the wavefunction overlap of two Majorana states along 
the junction is vanishingly small. For this reason, we 
are justified to treat these Majorana states as essentially 
non-interacting, except when a pair of them are near the 
edge of the TI. 

Near each 7r-crossing the zero-energy Majorana state is 
accompanied by other non-zero energy Andreev-bound- 
states (ABS's). Approximating the mass term by a linear 
potential, E(y) = i:A(y — i/i)/£b, (valid in the vicinity of 
9(y) « 7r), these ABS's have energies ±E n , where E n ss 
^/2-kuvmA/Ib, n = 1, 2, . . . , n max , and n max w £ B A/v M . 

Contrast to Conventional Junctions - Low-energy 
Andreev bound-states can also occur near ir phase dif- 
ference in conventional metallic Josephson junctions un- 
der ideal conditions. For a transparent superconductor- 
normal interface, and in the absence of spin-orbit cou- 
pling, a conventional 2D metallic junction can be thought 
of as two separate copies of Eq. 3 one for right mov- 
ing spin-up electrons and left-moving spin-down holes, 
and another for the corresponding time-reversed part- 
ners. This would lead to low-energy ABS's bound to 
Josephson vortices similar to those described above, ex- 
cept doubly degenerate. 

However, in practice normal scattering at the S-N in- 
terface (for example due to the chemical potential mis- 
match in S and N region) and spin-orbit coupling will 
push ABS up in energy towards the bulk gap A, and 
even completely remove them in a short Josephson junc- 
tion. In contrast, the helical nature of the TI surface 
state gives rise to topologically protected ABS's, which 
is also robust against disorder. Therefore, in the short- 
junction limit, the Josephson signatures described below 
are particular to the special properties of the S-TI-S junc- 
tion. 

Josephson Current - Having discussed the structure 
of low-energy Andreev bound-states in the junction, we 
now analyze their effect on Josephson current. When the 
number of magnetic flux quanta, $b/&0j is not close to 
a non-zero integer then the Josephson current is carried 
predominately by conventional states (extended states 
with E > A, and non-zero energy ABS's) and shows a 
nearly sinusoidal current phase dependence. Near integer 
values of flux quanta, &b/&o = ±l,±2,etc. . . , however, 
the conventional contribution becomes vanishingly small. 



In this regime, the Majorana bound-state contribu- 
tion to Josephson current dominates. For simplicity, we 
will first describe the situation for = $ in detail. 
The case of $ B = A$ is similar. When $ B = $ > 
and #o 7^ 7T, there is exactly one Josephson vortex 
piercing the junction, which binds Majorana modes at 
yo = ^f(l — ^j 1 ) on both the top and bottom surfaces 
of the TI. These Majorana modes hybridize by tunnel- 
ing into each other around the perimeter of the junction. 
The resulting energy splitting is exponentially suppressed 
as e -2min[y ,(w-3/o)]/AB 5 and is ne gligibly small when the 

position of Majorana modes yo is more than a few As 
away from the edges of the junction. In this regime, the 
splitting is insensitive to yo and 6q, and the Majorana 
modes do not contribute to the Josephson current. 

In contrast, when 9q w tt, yo w and the Majorana 
states are strongly coupled near the edge and split away 
from zero-energy. If the junction height h< \b, then we 
may ignore then finite thickness of the TI film, and the 

Majorana states are split by energy SEm ~ \J V i'^ ■ As 
0q approaches ir from below, the Majorana states move 
towards the junction edge at y — and begin to fuse 
and split when #o~7r(l — ^r). Increasing 9o beyond tt 
causes a different set of Majorana states to emerge from 
near y = W, and move to decreasing y, reversing the pro- 
cess that occurred near y = for 9 a > tt. The hybridiza- 
tion of Majorana states at the two edges gives rise to local 
supercurrents in opposite directions. At 9 = tt, the split- 
ting of Majorana states is large, but supercurrents from 
two edges cancel. Slight deviation from 9q — tt tips the 
balance by increasing the magnitude of supcrcurrent at 
one edge and suppressing the other, thereby generating a 
nonzero total supercurrent. The sensitivity of the Majo- 
rana splitting energy to the phase difference implies that 
the Majorana states contribute a peak in the Josephson 
current near 4> a w tt (l — ^) followed by a dip in the 
Josephson current near (f> w tt (l + 

By these considerations, we find that the maximal 
supercurrent carried by the Majorana state is Im ~ 
f^^f « A/$ , which for Al SC layers is I M ~ 10nA. 
We note that this result is largely independent of details. 
In particular it is completely independent of %, /x, W, 
and L (though it can depend slightly on finer details such 
as the degree of asymmetry between the top and bottom 
surfaces, or between the two edges). The Majorana cur- 
rent Im ~ A/<I>o roughly corresponds to the maximal 
amount of Josephson current carried by a single quan- 
tum channel[27]. By comparison, for junctions of a few 
am in width, such as those measured in Ref. 20, there are 
roughly kpW w I0 2 — I0 3 quantum channels in the en- 
tire junction. Consequently, the Majorana contribution 
to Josephson current, which dominates near an integer 
flux quanta, is 10~ 2 — 10~ 3 smaller than the maximum 
supercurrent for zero flux quanta. 

The situation is similar when a larger integer number, 
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FIG. 2. Numerical computation of the Andreev bound-state 
spectrum, and Josephson current for Eq. 3 with the param- 
eters v F = 4.2 x 10 5 m/s, fj, = lOmeV, A = 151//eV, and 
W = 2/im relevant for Al/Bi2Se3/Al junctions. Panels a and 
b show the low-lying Andreev bound-state spectrum as a func- 
tion of phase difference 9 between the superconductors for 
<E>s = 0.5<3?o and <E>s = <3>o respectively. The Josephson cur- 
rent corresponding to a) and b) are shown in panels c) and d) 
respectively. Panel d) displays the characteristic sharp peak- 
dip structure from topological Andreev bound-states fusing 
at the edge of the junction, as discussed in the text. 

N , of magnetic-flux quanta pierce the junction. Here, 
there are N Josephson vortices, each with a bound Ma- 
jorana state. For N <C j^, the splitting of these states 
from tunneling between Majoranas is negligible, except 
in the vicinity of 8q ~ ±7r where Majorana states fuse 
and annihilate at the edges of the junction. A similar 
peak-dip current-structure is observed for higher integer 
number of flux, N, with maximal current Im ~ 5^ inde- 
pendent of N, and with the peak (dip) width 56 ~ . 

Like the Majorana zero-modes, as t/i — > 0, W, the non- 
zero energy Andreev-bound-states on the top and bottom 
surfaces hybridize and split. However, unlike the Majo- 
rana modes, both branches of the hybridized finite-energy 
ABS's arc occupied and tend to give opposite and can- 
celling contributions to the Josephson current. 

Numerical Solution - To illustrate the predictions 
of the semiclassical treatment of this simple model, we 
have also numerically computed the spectrum of Andreev 
bound-states and their contribution to Josephson current 
by discretizing Eq. 3 on a finite lattice of spacing a, with 

E(y) w A cos {^p^j ■ In this simulations, we used realistic 

parameters of vf and A for an Al/Bi 2 Se3/Al junction, 
and junction geometry comparable to that in Ref. 20- 
22. While \x for Bi 2 Se3 is typically w O.leV, for simplic- 
ity, we illustrate the case for smaller fj, w lOmeV which 
can readily be achieved by gating[17, 22]. The Andreev 
bound-state spectrum and corresponding Josephson cur- 
rent are shown in Fig. 2 for half a magnetic flux quantum 
(left panel) and a single magnetic flux quantum (right). 
For half a magnetic flux quantum, there is a single 
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FIG. 3. Current-phase relationship for the parameters listed 
in Fig. 2 for a wide range of flux, (left) and for $b within 
a few percent of "l?o (right). Each curve is labeled by value of 
$s/^o- The Majorana contribution to the Josephson current 
becomes appear as a shoulder in the curves of the right panel 
which grows and eventually dominates very close to $ = $o- 
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FIG. 4. Critical current, Jo, as a function of flux, $b in linear 
(left) and log- scale (right). The curve is quite close to the 
conventional Frauenhofer pattern, with the exception that Ic 
does not vanish at integer number of flux quanta due to the 
extra contribution of the Majorana bound-states (as can be 
seen in the log-scale plot). 

Majorana zero- mode when 9q € [0,7r], however, the cur- 
rent is dominated by conventional states. For a full flux 
quantum, there are Majorana bound-states for all values 
of #0; which fuse at the junction edge when 6 w ir as 
described above. Away from 9 w ir, the spectrum is 
roughly independent of phase. Near 9 = ir, the ABS's 
near the 7r-crossings on the top and bottom surfaces hy- 
bridize and contribute to Josephson current. However, 
since the non-zero energy ABS's split in opposite direc- 
tions the Josephson current is predominately carried by 
the lowest-energy state, as described above. 

The results are similar for larger chemical potential, ex- 
cept that Vm is smaller, leading to a smaller level spacing 
for the non-zero energy ABS's, and a narrower peak-dip 
structure at $b = $0. Importantly though, the maxi- 
mum current in the peak-dip structure shown in Fig. 2d. 
is independent of vm , in agreement with the analytic ar- 
guments presented above. 

Fig. 3 shows the current phase-relation for various val- 
ues of flux through the junction. The Majorana mode 
contribution only starts to become visible for flux within 
s=s 5% of a single flux quantum (see right panel of Fig. 3). 
This contribution appears initially as a shoulder on the 
background sine-curve, which strengthens and becomes 
dominant within 1% of a single flux quantum. 

Finally, Fig. 4 shows the critical current, Ic, as a func- 
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FIG. 5. Current-phase relationship for a single magnetic flux, 
and varying TI thickness h = Onm, 20nm, 40nm, and 60nm. 
Thicker sides tend to contribute conventional Josephson sig- 
natures that mask the topological Andreev-bound-state con- 
tributions. 

tion of flux through the junction. Iq nearly follows the 
characteristic Frauenhofer pattern, except that Ic docs 
not vanish for integer flux due to the Majorana contri- 
butions to the Josephson current. In the experiments of 
Ref. 20, the avoided zero in the Frauenhofer pattern is 
substantially larger (w 10% of the maximal Ic)- As de- 
scribed below (and illustrated in Fig. 5), this discrepancy 
can be explained by accounting for the non-negligible 
thickness of the TI film in these experiments. 

Discussion - Up to now we have assumed that the 
thickness, h, of the TI film is smaller than the size of the 
Majorana states, A#. In practice, As is at most a few 
tens of nanometers. For thicker films with h > X B , the 
side of the TI may host additional states compared to 
h = when the phase difference at the edge is close to 
7r. As shown in Fig. 5, these states will conduct Joseph- 
son current, as their energies are sensitive to the phase 
difference along the sides of the junction. Like the Ma- 
jorana contribution described above, the contribution of 
current from these sides can dominate when there are 
close to an integer number of magnetic-flux quanta pierc- 
ing the junction. Unlike the Majorana contribution how- 
ever, these states exhibit a more conventional sinusoidal 
current-phase relationship, rather than a sharp peak-dip 
structure (see Fig. 5). 

Therefore, for the purpose of observing Josephson- 
current signatures of Majorana fermions, it is advanta- 
geous to make the TI film as thin as possible without 
strongly hybridizing the top and bottom surfaces. Since 
the confinement length in the z-direction is only a few 
nanometers, this last constraint is not too severe. Fur- 
thermore, it is advantageous to limit the junction width 
W (while maintaining W > Xb) in order to suppress the 
contributions from conventional extended states (which 
scales as I ~ W, in contrast to the width independent 
contribution from the Majorana bound-states). 

Lastly, while we have mainly considered the geometry 
shown in Fig. 1, where superconductors are deposited 
on the top and bottom surfaces. In practice it is easier 
to fabricate devices where superconductors are deposited 
only on the top surface. As described above, we expect 



our analysis to also apply to this simpler to fabricate ge- 
ometry, so long as the TI thickness is less than a coher- 
ence length, and so long as superconductivity is transmit- 
ted from top to bottom surfaces via either bulk TI states 
or boundary states on side surfaces. The main difference 
from the case discussed in the text is that there will be 
a non-zero hybridization of the ABS's on the top and 
bottom surfaces through the bulk. Nonetheless, in this 
case the energy splitting of the top and bottom ABS's is 
position independent, which will not qualitatively alter 
the Josephson signatures discussed above so long as the 
splitting is much smaller than the bulk gap A. 

Throughout this work we have implicitly assumed that 
the system is in its ground state. Finite temperature will 
not qualitatively alter our results, provided T < 6 Em- 
Deviations from ground-state behavior may also occur 
due to sources of single electrons, such as localized impu- 
rity states near the junction. When the Majorana states 
are far from the junction edge, their mutual fermion par- 
ity can be switched by tunneling to local single-electron 
sources. After such a parity switching, the Majorana 
modes will follow the positive energy branch as they 
approach the junction boundary and fuse, thereby con- 
tributing the opposite sign of Josephson current com- 
pared to the equilibrium case discussed above. Such 
parity switching events can thereby lead to hysteretic 
current-phase behavior, whose observation would provide 
strong evidence for the Majorana character of the An- 
dreev bound-states in the junction. 
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